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ABSTRACT 

The traditional first-order Comp ton-Get ting effect, which relates 
particle distributions as observed in two frames of reference moving 
with constant relative velocity, is inadequate for the description of 
low energy particles (less than a few hundred keV/nucleon) in the solar 
system. An exact procedure is given for recovering both isotropic and 
anisotropic distributions in the solar wind frame from observations 
made in a spacecraft frame. The method is illustrated by analyzing a 
particle event observed by the University of Maryland experiment on 


IMP-7 on 31 October 1972. 
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INTRODUCTION 

A distribution of particles which is isotropic in one frame 
of reference will display an anisotropy if observed from a different 
frame of reference. This is referred to as the Comp ton-Get ting effect 
(Compton and Getting, 1935). For nonrelativistic particles the magni- 
tude of the induced anisotropy (Gleeson and Axford, 1968; Forman, 1970) 
is equal to (2 + 2y) (w/v) , where y is the spectral index of the parti- 
cle differential intensity, w is the relative speed of the two frames 
of reference, and v is the particle speed. The above formula assumes 
w << V. Recently Balogh, et al. (1973) have derived an expression for 
the anisotropy accurate to order (w/v)^. For particle convection induced 
by the solar wind, w is equal to the solar wind speed. Present day 
satellite experiments respond to such low energies that the assumption 
w << v can no longer be used. For example, one detector on the Univer- 
sity of Maryland experiment on IMP-8 responds to heavy ions with energies 
v20 keV/nucleon (Turns, et al. , 1974). Since v (km/sec) = 440 [E(keV/ 
nucleon) ] ^ , this implies v 5 w for typical solar wind speeds , and 
V = 3 w during disturbed periods. Future experiments will undoubtedly 
reach even lower energies. 

In this letter we first derive an exact transformation of a particle 
distribution which is isotropic in the solar wind frame. We then show 
how a distribution with arbitrarily high anisotropy in the solar wind 
frame may be transformed exactly into the observer's frame. The procedure 
is illustrated by analyzing one of the "post-shock" particle spikes dis- 
cussed by Gloeckler, et al (1974) . 



EXACT TRANSFORMATION FOR ISOTROPTC DISTRIBUTION 


The exact transformation procedure is based on the Lorentz invari- 
ance of the particle distribution fvmction in phase space (Forinan, 1970). 
Let primed quantities refer to ts.c solar wind frame and unprimed quan- 
tities to the observer's frame. All particles arc assumed to be non- 
relativistlc, since for relativistic particles the first-order Compton- 
Getting correction is a perfectly good approximation* 

The particle momenta in the two frames of reference are related by 
the Lorentz transformation 

P * = P “ P w/v * (1) 

Here w is the solar wind velocity and v is the particle speed. (Equa- 
tion (1) assumes w << speed of light) . The magnitudes of the momenta 
are related by 

p’ = P [1-2 (w/v) cos 0 -f (w/v)^]^/^, (2) 

where 0 is the angle between the solar wind velocity and the direction 

.-V 

in which the observer is looking (i.e,, cos 9 = (w*P)/(wP)). 

The distribution function f(P) is the nurnber of particles of momen- 
tum ? in the volume element dP^, ‘iPyj phase space. 

As demonstrated by Forman (1970) , this function is a Lorentz invariant: 



4 


For the case of a distribution which is isotropic in the solar wind 
frame we have 

f(P) = f ' (P’) . (4) 


If in addition the distribution fxmction is a power law in the solar 
wind frame , f'(P’) “ (P')~ , then equations (2) and (4) may be, com- 
bined to yield 

f(P) = f’(P) [l-2(w/v) cos 0 + (w/v)^] . (5) 


It is apparent from equation (5) that the spectrum is not a power law . 
in the observer's frame; both the shape and the intensity of the spec- 
trum are functions of angle and energy. This may be seen quantitatively 
by computing an effective spectral index for a given momentum and angle 
as : 




- <i(ln f) 

" " d(ln P) 


( 6 ) ' 


Using equation (5) this becomes 


„ = „t r 1 - (w/v) cos Q , 

%ff ^ 1-2 (w/v) cos 0 + 6^/v) ^ 


(7) 


Introducing the particle differential intensity dJ/dE (particles per 
sec - ster ^ keV) and using the relation dJ/dE = P^f (Forman, 1970), 
it follows that if the distribution function is a power law in momentum, 
f « P then the differential intensity is a power law in kinetic energy, 
dJ/dE cc with y = (n/2) - 1. Thus, using equation (7), the effective 
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differential intensity spectral index is 

, (y’+1)(w/v) [cos 0 - (w/v)] ... 

^eff ~ ^ 1-2 (w/v) cos 0 + (w/v)'^ 

Figure 1 shows plotted vs. 0 for (v/w) = 4 and 10, assuming a power 

law spectrimi In the solar wind frame with y* “ 3. The steepest spectrum 
is seen when looking towards the sun (0 = 0'^). 

Again using the relation dJ/dE = obtain, using equation (5), 

= const X E ^ [1-2 (w/v) cos 0 + (w/v)^3 ^ • (9) 

qE 

This differential intensity is plotted in Figure 2 for three values of 0 
assuming a solar wind speed of w = 400 km/sec and a spectral index - 3 
in the solar wind frame. 

For small (w/v) we can expand equation (5) using n = 2y 4 2 to 
obtain, correct to order (w/v)^, 

f(?) == f’(p) [1 + C (w/v) + D (w/v) 2] , (10) 


where 


C = 2(y’ + 1) cos 0 

D = (y' + 1)2 + (y’ + 1) (y* + 2) cos 20 

Here C is the familiar first-order Compton-Getting correction. The second 
order term, D, has been obtained previously by Balogh, et al. (1973) by 
expanding equation (4) in a Taylor series. 



6 


EXACT TRANSFORMATION FOR ANISOTROPIC DISTRIBUTION 

We describe both frames of reference using spherical coordinate 
systems with the polar axes along the sun earth line, and again use 
primes to denote quantities defined in the solar wind frame. Thus the 
momentam P in the observer's frame is labeled by magnitude P, polar 
angle 0, and azimuthal angle <ji. Since the azimuthal angles span planes 
which are normal to the solar wind velocity, it follows that = 4>’ 

(i.e. , in cartesian coordinates P^, = » ^y ~ ^y'^’ change in the 

magnitude of momentum is given by equation (2) . The polar angle G' must 
also be Lo rents- transformed (see e.g., Jackson, 1962): 


„ _i r sin 0 1 

® = Tan M cos 0 - Tw'/v)" ^ 


( 11 ) 


Thus equation (3) can be written quite generally as 


f(P,0,it)) = f ' (P[l-2(w/v)cos0 + (w/v) 2]1/2^ Tan'^ ^ cosl ^ ^^vf ^ 


In order to illustrate the use of equation (12) , it is convenient to 
make a number of simplifying assumptions* Firstly, we assume that the 
detector system scans in the ecliptic plane so that the angle ^ need not 
be considered* (It is a straightforward matter to consider different 
orientations.) Secondly, we assume that the distribution function in the 
solar wind frame is separable in the sense that 


f'(P’) = g'(P') h’(0’) * 


(13) 
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This is probably a reasonable assumption over a limited energy range, 
although it should be noted that equation (2) implies even a detector 
responding only to an energy will sample particle energies in the 
solar wind frame ranging from [1 - (w/v)]^E^ to [1 + (w/v)]^E^. Of 
course, if equation (13) does not allow a sufficiently good fit to the 
experimental data, a more general form must be used. Finally, we 
assume a power law spectrum with g’(P’) “ (P’) ^ . Now equations (5), 
(12), and (13) may be combined to yield the differential intensity seen 
in the observer’s frame: 


dE 


cons 


t X E~^ [1-2 (w/v) cos 6 + (w/v)2] 


xh’ {Tan-1 [ • 

*^cos 0 - (w/v) 


sin 0 


(14) 


APPLICATION 

The use of equation (14) will be illustrated by analyzing one of 
the so-called "post-shock" spikes observed in interplanetary space by the 
University of Maryland experiment on IMP-7, and discussed by Gloeckler, 
et al. (1974) and Levy, et al. (1974). This particular sharp intensity 
increase lasted from 1954 to 1956 UT on October 31, 1972. The detector 
we are concerned with responds to alpha particles with energies between 
95 and 135 keV/nucleon. The rate is sectored into four 90° sectors in the 
ecliptic plane, so that it is not possible to uniquely recover the func- 
tion h' in equation (14). (This indicates the importance of additional 
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sectoring for future low energy experiments.) Instead, we must arbitrar- 
ily assume a form for h’, numerically integrate equation (14) over four 
90® sectors, and then compare the resultant four numbers with the four 
observed rates. Both observed and calculated rates are normalized to give 
an average rate of 50 for each sector; there are thus three independent 
rates to fit. We have assumed a two-parameter particle distribution in 

the solar wind frame of the form h'(0') = 1 + b cos (O' - 0 ') with b 

o ’ 

representing the magnitude of the anisotropy and 0^’ its direction. 

The solar wind velocity during this time was w ~ 700 km/sec (H. 
Rosenbauer and H. Grunwald, private communication, 1974), so that v = 7w. 
The average differential intensity spectral index, obtained from the 
observed intensities summed over all sectors for two detectors at differ- 
ent energies, was y = 3.0 from 'viOO to 400 keV/nucleon. From Figure 1, 
we expect this average value to be approximately equal to the value in 
the solar wind frame, so we take y’ = 3.0. Thus the normalized rates 
obtained by integrating equation (14) depend on only two variables; b and 
Therefore, if the three independent observed rates can be fit, some 
confidence can be placed in the assumed form for h*. In this example, 
the best fit was obtained by choosing b = 0.95 and 0^' = 146°. In Figure 
3 the dotted curve represents the assumed angular distribution of particles 
in the solar wind frame. The solid curve (obtained from equation (14)) 
represents the same distribution as viewed from the observer’s frame. The 
normalized observed rate (with one sigma error) and the calculated rate 
are shown for each sector. 
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Since the agreement is really quite good, we conclude that our 

assumed form for h' is an accurate representation of the actual particle 

distribution in the solar wind frame during this event* Although the 

magnitude of the derived anisotropy is large (95%) , the particle distri'** 

bution is broad - 177^ in the solar wind frame)* A more sharply 

peaked distribution cannot fit the data* This may be seen qualitatively 

by considering a delta function distribution in the solar wind frame: 

h’(9’) =6(0’ - 0 ’). Equation (14) shows that the distribution would 
o 

also be a delta function in the observer’s frame, 6(0 - * Equation 

(11) shows that for 0 ’ - 146°, G - 141°. Thus the calculated rate 
would be 200 for sector 1, and zero for the other three sectors. This 
would obviously be a poor fit to the observed rates shown in Figure 3, 

The average direction of the inagnetic field in the ecliptic plane 
during this time period, observed on the same satellite, was 152° ± 6° 
(shown by the arrow in Figure 3) , where the variation represents real 
fluctuations in the field direction (N* Ness and R. Lepping, private 
communication, 1974)* Thus the particle streaming was essentially field 
aligned with the particles coming from an easterly, anti -solar direction. 
From similar analyses, this same conclusion also holds for the other 
post-shock spikes described by Gloeckler, et al* (1974), 
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FIGURE CAPTIONS 


Figure 1 . Spectral index in the observer's frame as a function of 
direction assuming a spectral index y' = 3 in the solar wind frame. 

Figure 2. Differential intensity in the observer's frame for three dif- 
ferent directions assuming a solar wind speed of 400 km/sec and an iso- 
tropic power law differential intensity in the solar wind frame with 
= 3. 

Figure 3 . Angular particle distributions in the solar wind frame (dotted 
curve) and satellite frame (solid curve) deduced from observations of 
alpha particles (95-135 keV/nucleon) in a ’’post-shock" spike (1954-1956 
UT October 31, 1972)* The upper number in each bracket is the observed 
normalized counting rate (with one sigma error) for each sector, and the 
lower number is the calculated rate obtained by integrating the solid 
curve over that sector. The arrow indicates the average direction of the 
magnetic field in the ecliptic plane during this time period. 




Figure 1 





